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We study the cosmological and weak-field properties of theories of gravity derived by extending 
general relativity by means of a Lagrangian proportional to R 1+s . This scale- free extension reduces 
to general relativity when 8 — » 0. In order to constrain generalisations of general relativity of 
this power class we analyse the behaviour of the perfect-fluid Friedmann universes and isolate the 
physically relevant models of zero curvature. A stable matter-dominated period of evolution requires 
8 > or 8 < —1/4. The stable attractors of the evolution are found. By considering the synthesis 
of light elements (helium-4, deuterium and lithium-7) we obtain the bound —0.017 < 8 < 0.0012. 
We evaluate the effect on the power spectrum of clustering via the shift in the epoch of matter- 
radiation equality. The horizon size at matter-radiation equality will be shifted by ~ 1% for a value 
of 8 ~ 0.0005. We study the stable extensions of the Schwarzschild solution in these theories and 
calculate the timelike and null geodesies. No significant bounds arise from null geodesic effects but 
the perihelion precession observations lead to the strong bound 8 — 2.7 ±4.5 x 10~ 19 assuming that 
Mercury follows a timelike geodesic. The combination of these observational constraints leads to 
the overall bound < 8 < 7.2 x 10~ 19 on theories of this type. 

PACS numbers: 95.30.Sf, 98.80.Jk, 04.80.Cc, 98.80.Bp, 98.80.Ft, 95.10.Eg 



I. INTRODUCTION 



There is a long history of considering generalisations of Einstein's theory of general relativity which reduce to 
general relativity in the weak gravity limit when the spacetime curvature, i?, becomes small. Typically, these studies 
consider a gravitational Lagrangian which augments the linear Einstein- Hilbert Lagrangian by the addition of terms 
of quadratic or higher order in R, first considered by Eddington p| ; these additions may also include terms in In R, 
. More general extensions of general relativity in this spirit have considered the structure of theories derived from 
gravitational Lagrangians that are general analytic functions of R, |j.Lll0@ These choices produce theories which 
can look like general relativity plus small polynomial corrections in the appropriate limiting situations as R becomes 
small. There has also been interest in theories with corrections to general relativity that are 0(R~ l ) because of their 
scope to introduce cosmological deviations from general relativity at late times which might mimic the effects of dark 
energy on the Hubble flow 0, Q . We also know that theories derived from a Lagrangian that is an analytic function 
of R have an important conformal relationship to general relativity with scalar field sources so long as the trace of the 
energy-momentum tensor vanishes in the higher-order gravity theory 0, ll(| . All these theories introduce corrections 
to general relativity which come with a characteristic length scale that is determined by the new coupling constant 
that couples the higher-order terms to the Einstein-Hilbert part of the Lagrangian. In general, these theories are 
mathematically complicated with 4 t,l -order field equations that can exhibit singular perturbation behaviour unless 
care is taken to ensure that the stationary action does not become maximal rather than minimal 0, 1 1 1| . and there 
are few interesting exact solutions other than those of general relativity, which are particular solutions in vacuum and 
for trace-free energy momentum tensor so long as the cosmological constant vanishes [f| . 

In this paper we are going to consider a different type of generalisation of Einstein's general relativity, in which no 
new scale is introduced. The Lagrangian is proportional to R n , and so general relativity is recovered in the n — ► 1 
limit, from above or below. Particular cases have been studied by Buchdahl an d Roxburgh [l^. This gravitation 
theory has many appealing properties and, unlike other higher-order gravity theories, admits simple exact solutions for 
Friedmann cosmological models and exact static spherically symmetric solutions which generalise the Schwarzschild 
metric. As well as allowing comparison with observation these solutions also provide an interesting testing ground for 
new developments in gravitation theory such as particle production, holography and gravitational thermodynamics. 
Furthermore, this theory is of additional interest because it permits a very general investigation of the nature of its 
behaviour in the vicinity of a cosmological singularity which brings the behaviour of general relativity into sharper 
focus. In another paper we show that the counterpart of the Kasner anisotropic vacuum cosmology can be found 
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exactly and strong conclusions drawn about the presence or absence of the chaotic behaviour found in the Mixmaster 
universe. 

The structure of this paper is as follows; in the next section we present the gravitational action and field equations for 
the theory of gravity we will be considering. A conformal relationship with general relativity containing a scalar field 
in a Liouville (exponential) potential is then outlined and the Newtonian limit of the field equations is investigated. 
The rest of the paper is then split into two sections; the first investigates the cosmology of the theory and the second 
investigates the static and spherically symmetric weak field - in both cases our goal is to calculate predictions for 
physical processes, the results of which can be compared with observation. We use observational data from cosmology 
and the standard solar system tests of general relativity to bound the allowed values of n, the single defining parameter 
of the theory. 

In the cosmology section we consider Friedmann-Robertson- Walker universes. We present the equivalent of the 
Fricdmann equations, in this theory, and find some power-law exact solutions. A dynamical systems approach is 
then used to show the extent to which these solutions can be considered as attractors of spatially flat universes at 
late times. After showing the attractive properties of these solutions (with certain exceptions) we proceed to predict 
the results of primordial nucleosynthesis and the form of the power spectrum of perturbations in this theory. These 
predictions are then compared to observation and used to constrain deviations from general relativity. 

The static and spherically symmetric weak-field analysis follows. We present the field equations and find the 
physically relevant exact solution to them. A dynamical systems approach is then used to find the asymptotic 
attractor of the general solution at large distances. This asymptotic form is then perturbed and the linearised field 
equations are found and solved. The exact solution is shown to be the relevant solution in this limit, when oscillatory 
modes in the perturbed metric functions are set to zero. We find the null and time-like geodesies for this spacetime 
to Newtonian and post-Newtonian order. Predictions are then made for the outcomes of the classical tests of general 
relativity in this theory; namely the bending of light, the time-delay of radio signals and the perihelion precession of 
Mercury. These predictions are then compared to observation and again used to constrain deviations from general 
relativity. 

II. FIELD EQUATIONS 

We consider here a gravitational theory derived from the Lagrangian density 

1 —,1+5 



C G = -V~gR l+ \ (i) 

x 

where 6 is a real number and x is a constant. The limit 5 — *■ gives us the familiar Einstein-Hilbert Lagrangian of 
general relativity and we are interested in the observational consequences of |<5| > 0. 
We denote the matter action as S m and ignore the boundary term. Extremizing 

Ccd 4 x + S m , 



with respect to the metric g ab then gives [lj] 

S(i S 2 )R S ^± - 8(1 + 8)R S ^ + (1 + 8)R s R ab - \g ab RR 5 



- g ab S(l - 8 2 )R A —^- + 8(1 + 8)g ab R A — = ^-T ab , (2) 
R A R 2 

where T ab is the energy-momentum tensor of the matter, and is defined in terms of S m and g ab in the usual way. We 
take the quantity R s to be the positive real root of R throughout this paper. 

A. Conformal equivalence to general relativity 

Rescaling the metric by the conformal factor O(r) = QqR s the vacuum field equations become 

= 3d 2 R. a R, b 3<5 2 _ _ cd R }C R,d 8 g a b R 

Lrofe — -T, 7^, T9abg 



2 R 2 4 yaDJ R 2 2(1 + 8) n R S1 

where g ab = £lg ab and other quantities with overbars are constructed from the rescaled metric g ab . 
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Making the definition of a scalar field 



:! lni? 5 , 



16ttG 



these equations can be rewritten as 



and 



where V(<f>) is given by 



a,,, ■■ s, T (V ( <f> >a 4>, b - \ga»{9 cd <i>A,d + 2V&)) ) 0!) 



U( P = "IT' 



T _... <5sign(i?) /16ttG(1-5), 

^ = 16^(1 + ^ eXP \V- T 1 *) ■ (4) 

The magnitude of the quantity Slo is not physically important and simply corresponds to the rescaling of the metric 
by a constant quantity, which can be absorbed by an appropriate rescaling of units. It is, however, important to 
ensure that Sl > in order to maintain the +2 signature of the metric. This result is a particular example of the 
general conformal equivalence to general relativity plus a scalar field for Lagrangians of the form f(R), where / is an 
analytic function found in refs. 0,0- 

B. The Newtonian Limit 

By comparing the geodesic equation to Newton's gravitational force law it can be seen that, as usual, 

r&, - *,m (5) 

where $ is the Newtonian gravitational potential. All the other Christoffel symbols have T a bc — 0, to the required 
order of accuracy. 

We now seek an approximation to the field equations that is of the form of Poisson's equation; this will allow us 
to fix the constant \- Constructing the components of the Riemann tensor from JSJ we obtain the standard results 

The 00 component of the field equations J2J can now be written 

(l + Wo-^ooi?=f^ (7) 

where terms containing derivatives of R have been discarded as they will contain third and fourth derivatives of $, 
which will have no counterparts in Poisson's equation. Subtracting the trace of equation Q gives 

( 1 + Wo = ^(Too-^ y5 oor) (8) 

where T is the trace of the stress-energy tensor. Assuming a perfect-fluid form for T we should have, to first-order, 

T 00 - P and T ~ 3p - p ~ -p. (9) 

Substituting © and © into JSJ gives 



4(1 -P) R s ' 
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Comparison of this expression with Poisson's equation allows one to read off 

X = 16ttG jj— Qr 5 (10) 

where Rq is the value of the Ricci tensor at the time G is measured. It can be seen that the Newtonian limit of the 
field equations J2J does not reduce to the usual relation V 2< £ oc p, but instead contains an extra factor of R 5 . This can 
be interpreted as being the space-time dependence of Newton's constant, in this theory. Such a dependence should be 
expected as the Lagrangian (0 can be shown to be equivalent to a scalar-tensor theory, after an appropriate Legendre 
transformation [46j (see e.g, |16|). This type of Newtonian gravity theory admits a range of simple exact solutions in 
the case where the effective value of G is a power-law in time even though the theory is non-conservative and there 
is no longer an energy integral [Tsj . 



III. COSMOLOGY 



In this paper we will be concerned with the idealised homogeneous and isotropic space-times described by the 
Friedmann-Robertson-Walker metric with curvature parameter k: 

ds 2 = -dt 2 + a 2 (t) ( - ^ 9N + r 2 d9 2 + r 2 sin 2 OdcjA . (11) 
\ (1 — KT 1 ) J 

Substituting this metric ansatz into the field equations J5J , and assuming the universe to be filled with a perfect fluid 
of pressure p and density p, gives the generalised version of the Friedmann equations 

(1 - 5)R'+ S + 35(1 + S)R S (J + - 35(1 - 6 2 )R S ^ = |(p - 3p) (12) 

-3-(l + S)R S + ^ + 35(1 + 5)~R S = f p (13) 
a I a R 2 



where, as usual, 



i? = 6-+6^+64. (14) 



It can be seen that in the limit (5^0 these equations reduce to the standard Friedmann equations of general relativity. 
A study of the vacuum solutions to these equations for all k has been made by Schmidt, see the review |17| and a 
qualitative study of the perfect-fluid evolution for all n has been made by Carloni et al Various conclusions 

are also immediate from the general analysis of f(R) Lagrangians made in ref p| by specialising them to the case 
/ = R 1+s . In what follows we shall be interested in extracting the physically relevant aspects of the general evolution 
so that observational bounds can be placed on the allowed values of 6. 

Assuming a perfect-fluid equation of state of the form p = cap gives the usual conservation equation p cx a^ 3 ^ +1 \ 
Substituting this into equations (|12|l and l|13fl . with k = 0, gives the power-law exact Friedmann solution for 

2(1 + 6) 

a(t)=t^+nn (15) 

where 

(1 - 2S)(2 - 3(5(1 + w) - 25 2 (4 + 3w)) = 12ttG(1 -5)(1 + wfp c (16) 

and p c is the critical density of the universe. 

Alternatively, if u> = —1, there exists the de Sitter solution 



a(t) 



e 



where 



3(1 - 25)n 2 = 8vrG(l - 5)p c . 
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FIG. 1: Critical density, Qq, as a function of 5. Solid line corresponds to pressure-less dust and dashed line to black-body 
radiation. 



The critical density 1161) is shown graphically, in figure in terms of the density parameter fio = 8 oS c as a 

function of S for pressureless dust (u> = 0) and black-body radiation (u> = 1/3). It can be seen from the graph that 
the density of matter required for a flat universe is dramatically reduced for positive S, or large negative S. In order 
for the critical density to correspond to a positive matter density we require S to lie in the range 

V73 + 66w + 9w 2 + 3(1 + lu) . V73 + 66lu + 9uj 2 - 3(1 + u) 
< 5 < -. (17) 

4(4 + 3w) 4(4 + 3w) ^ ' 



A. The dynamical systems approach 

The system of equations ()12[) and 1|13|) have been studied previously using a dynamical systems approach by Carloni, 
Dunsby, Capozziello and Troisi for general k [l9|. We elaborate on their work by studying in detail the spatially flat, 
K = 0, subspace of solutions. This allows us to draw conclusions about the asymptotic solutions of l|12l) and 113|) 
when k — and so investigate the stability of the power-law exact solution 115|) and the extent to which it can be 
considered an attractor solution. By restricting to k = we avoid 'instabilities' associated with the curvature which 
are already present in general rclativistic cosmologies. 

In performing this analysis we choose to work in the conformal time coordinate 



*-Vs?* < 18 > 

Making the definitions 

R' , a' 

x — — and y = — , 
R y a ' 

where a prime indicates differentiation with respect to r, the field equations (|12f) and l|13|) can be written as the 
autonomous set of first order equations 

,_ 2-(5(l + 3^) Sx 2 (4-S(l + 3w))xy 2(1 - S)y 2 
X 5 2 {l + 5) 2 26 S 2 [ ' 

These coordinate definitions are closely related to those chosen by Holden and Wands [20[ for their phase-plane 
analysis of Brans-Dicke cosmologies and allow us to proceed in a similar fashion. 



1. Locating the critical points 

The critical points at finite distances in the system of equations H19fl and (|20fl are located at 

In = ± , and Ui 2 = ± — , (21) 

Sy/(l + S)(2-3w) ' ^/(l + 6)(2-3u) 
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and at 

T 3y/2(l+^) ygg + g 

134 =T — ; and w 3 4 = ± — . (22) 
* ^/(l + 5)(2 - 35(1 + <5 2 (8 + 6lo)) y ^J2 - 35(1 + to) - <5 2 (8 + 6w) 

The exact form of a(t) at these critical points, and the stability of these solutions, can be easily deduced. At the 
critical point (xi,yi) the forms of a(r) and -R(t) are given by 

o(r) = a e^ r and i?(r) = #oe XiT , (23) 

where ao and i?o are constants of integration. In terms of r the perfect-fluid conservation equation can be integrated 
to give 

where po is another positive constant. Substituting into the definition of r now gives 

dr oc e -U 1+u)v * T -i* tT db 

or, integrating, 

t-t oc^ - e i(.l+»)VtT+ix t r ( 24) 

Kl + wjyi + fii 

It can now be seen that if 3(1 + + fe, > then t — ► 00 as r — ► 00 and i -» ^ as r -» —00. Conversely, if 
3(1 + + 5xi < then t — > to as r — > 00 and t — > —00 as r — > —00. 

In terms of t time the solutions corresponding to the critical points at finite distances can now be written as 

a(t) oc (t - t ) 3 ( 1 +")»i+ i -i and R(t) oc (t - t Q )^+^n+^H . 

The critical points 1 and 2 can now been seen to correspond to a oc 1 2 and the points 3 and 4 correspond to l|15(l . 

In order to analyse the behaviour of the solutions as they approach infinity it is convenient to transform to the 
polar coordinates 

x = r cos 
y = f sin0. 

The infinite phase plane can then be compacted into a finite size by introducing the coordinate 

f 



1 + f 

The equations (|19|l and l|2(J|) then become 
-1 



45 2 (l + 5) 



4(1 - 2r)(£(l + 5) sin - (2 - 5(1 + 3w)) cos 0) 

- r 2 ((6 - 4(5 + 3<S 2 + 5 3 - 125w) cos0 + (1 + <5)(2 - 25 - 5 2 - 25 3 ) cos 30 

- 25(3 - 5(1 + 3uj) + 3 cos 20) sin 0)j (25) 

and 

(25(1 + 5) cos0 + 2(2 - 5(1 + 3u)) sin0)(l - 2r) 

- (5(1 + 5) cos 0(1-3 cos 20) - 4 sin + 4(1 - J) 2 sin 3 

+ 25(1 + 3u + 5(1 + 5)(1 + 25) cos 2 0) sin 0)r 2 . (26) 



25 2 (l + 5)(l-r)r 
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In the limit r — > 1 (f — > oo) it can be seen that critical points at infinity satisfy 

sin fa (8 cos <fii + sin<pi)(S(l + 28) cos<pi + 2(1 — 8) sm<fii) = 

and so are located at 

05,(6) = (+tt) (27) 
07,(8) =tan- 1 (-<5) (+tt) (28) 

09,(io)= tan- 1 (-f^) (+-)• (M) 

The form of a(t) can now be calculated for each of these critical points by proceeding as Holden and Wands [20| . 
Firstly, as r — » 1 equation l|25|l approaches 

r'-> 4^1 5(1 + 2(5(1 + 3cj)) sin 4 - 35 sin 3^ 
4d^ \ 

- (2 - 6(2 + S)) cos ^ + (2- 6(2 + 6+ 25 2 )) cos 3&J 

= f(4>i) 

which allows the integral 

r-l = /(^)(r-r ) 

where the constant of integration, tq, has been set so that r — > 1 as r — > To. Now the definition of x allows us to write 

R' r f((j)i) (t - t ) + 1 cos^j 

- ■ COS Qi — —, , ; ; ; — COS r 



i? (1-r) ri /(0i)(r-ro) ~ /(0i)(r-ro) 

as t — > To- Integrating this it can be seen that 

COB <fri 

ii oc It — tq| fWO as r — > 1. 



Similarly, 

The definition of t l|18|) now gives 
which integrates to 

where 



a oc |t — tq | /< * i) as J — ► 1. 



dT cx |t - t | 5( + w) 7<*7) + 27(*77^ 



t-iooc-^lT-TofW (30) 



_ 3(1 + lj) sin (pj +8 cos 0j - 2/(^) 

The location of critical points at infinity can now be written in terms of t as the power-law solutions 

R(t) oc (t-t )^*3 and a(t) oc (t - 1 )^^ . (31) 
Direct substitution of the critical points (|2Tjl . (|2~8l and l(2"5|) into l(3*Tl) gives 

®5,6(t) ~> constant 



a,7, s(t) — > Vt -t 

5(1 + 2,5) 

ag.ioW -» (*-*o) 

as r — > 1. Moreover, it can be seen from H3()|) that as r — ► 1 and t — > To so i — * £q as l° n g as F (fa) / f (fa) < 0) as i s 
the case for the stationary points considered here (as long as the value of (5 lies within the range given by l|17|l). 
The exact forms of a(t) at all the critical points are summarised in the table below. 
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Critical point 


a(t) 


1, 2, 7 and 8 


ti 


2(1+6) 


3 and 4 


£3(1+") 


5 and 6 


constant 




3(1 + 25) 


9 and 10 


t ci-«) 



,2. Stability of the critical points 

The stability of the critical points at finite distances can be established by perturbing x and y as 

x(r) — Xi + u(r) and y(r) = yi + v(r) 

and checking the sign of the eigenvalues, Xi, of the linearised equations 

u = XiU and v = \v. 

Substituting (1M2I) into equations (|19|l and l|20() and linearising in u and v gives 

/ (4-5(1 + 3c)) \ / (4-j(l + 3aQ) (1 + 5) \ 

u = - \ 6xi-\ — Vi )u - -r a;, + 4 — — — yi v 



(32) 



25 



, (2 + 35) / (2 + 3<5) 

«' = - — - — -2/0" + I - — - — -Xi 



25 

(2 + 5(1 + 3u)) 



5 2 



-2/o v 



2 " V 2 5 
The eigenvalues Ai are therefore the roots of the quadratic equation 

A 2 + BXi + C = 

where 

B = -i(2 + 5) Xi -|(l + 3w) W 

C = --(2 + 35):e. 2 - (2 + 5(1 + 35))^ + i(2 - 115 - 6^(1 - 5) + 95w 2 )yf. 
2 2d 

If B > and C > then both values of Xi are negative, and we have a stable critical point. If B < and C > both 
values of Xi are positive, and the critical point is unstable to perturbations. C > gives a saddle-point. 
For points 1 (upper branch) and 2 (lower branch) this gives 



B = T 



(I + 5(2 + 3lu) - 3lu) 



8y/(l + S)(2-3u) 
and for points 3 (upper branch) and 4 (lower branch) 

3(1-^(1 + 25)) 



and 



G 



B = ±- 



^2(1 + 5)(2 - 35(1 + u) - 25 2 (4 + 3w)) 



and 



(1 + 45 - 3oj) 
5(1 + 5) 

(l + 45-3w) 



C = 



5(1 + 5) 



The stability of the critical points at finite distances for a universe filled with pressureless dust are therefore, for 
various different values of 5, given by 



Critical point 


B 


c 


V73+3 < 5 < _ 


_i _i < 5 <0 0<5<^ 


1 


(1+28) 


(1+45) 


Saddle 


Stable 


Saddle 


5^2(1+5) 


.5(1+5) 


2 


(1+25) 


(1+45) 


Saddle 


Unstable 


Saddle 


V 2 (l+<5) 


5(1+5) 


3 


3 


(1+45) 


Stable 


Saddle 


Stable 


^2(l+5)(2-35-85 2 ) 


5(1+5) 


4 


3 


(1+45) 


Unstable 


Saddle 


Unstable 


\/2(l+5)(2-35-85 2 ) 


5(1+5) 



and for a universe filled with black-body radiation are given by 
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Critical point 



1 

2 
3 

4 




Ve+i ^ x ^ Ve-i 

_5 ^ _ 5 



Values of 6 < 



\Z73+66^+9^ 2 +3(l+oQ j 
4(4+3^) anC1 



A /(l+<5)(l-2<5-5(5 2 ) IT+^) 
V73+66a;+9a; 2 -3(l+u) 



s > 



Saddle 
Saddle 
Stable 

Unstable 

4(4+3^) have not been considered here as they lead to 

negative values of po for the solution (|15J) . (The reader may note the difference here between the range of 6 for which 
point 3 is a stable attractor compared with the analysis of Carloni et. al.). 

Point 3 lies in the y > region and so corresponds to the expanding power-law solution l|15p. It can be seen from 
the table above that this solution is stable for certain ranges of 6 and a saddle-point for others. In contrast, point 4, 
the contracting power-law solution (|15|) . is unstable or a saddle-point. The nature of the stability of these points and 
the trajectories which are attracted towards them will be explained further in the next section. 

A similar analysis can be performed for the critical points at infinity. This time only the variable (f> will be perturbed 

as 



fat) = fa+q(t). 



(33) 



The conditions for stability of the critical points are now that r' > and the eigenvalue p of the linearised equation 
q' = pq satisfies p < 0, in the limit r — > 1. If both of these conditions are satisfied then the point is a stable attractor, 
if only one is satisfied the point is a saddle-point and if neither arc satisfied then the point is repulsive. 
Substituting l|33|l into (|26[1 and linearising in q(t) gives, in the limit r — > 1, 



1 



q = A5 2 {l _ r) ((6(1 ~S) + 5 2 {l + 28)) cos fa - 3(2(1 - 5) - 8 2 {l + 26)) 



cos 3fa — 36 sin fa + 96 sin 3(f>) q 



The sign of r' in the limit r — > 1 can be read off from (|25|l . The stability properties of each of the stationary points 
at infinity can now be summarised in the table below 



Critical point 


-Ni <6<-± 


-i<<5<00<<S<i 


\ < 6 < N 2 


5 


Stable 


Saddle 


Unstable 


Unstable 


6 


Unstable 


Saddle 


Stable 


Stable 


7 


Unstable 


Unstable 


Saddle 


Stable 


8 


Stable 


Stable 


Saddle 


Unstable 


9 


Saddle 


Stable 


Stable 


Saddle 


10 


Saddle 


Unstable 


Unstable 


Saddle 



where N± 



_ V73+66^+9^ 2 +3(l+o;) 
~ 4(4+3w) 



and N 2 



\/73+66^+9a; 2 -3(l+cj) 
4(4+3w) 



3. Illustration of the phase plane 

Some representative illustrations of the phase plane are now presented. Firstly the compactified phase plane for a 
universe filled with pressureless dust, lo = 0, and a value of <5 = 0.1 is shown in figure|21 Figure|21is seen to be split 
into three separate regions labelled I, II and III. The boundaries between these regions are the sub-manifolds R = 0. 
As pointed out by Carloni et. al. the plane R — is an invariant sub-manifold of the phase space through which 
trajectories cannot pass. 

The equation for R in a FRW universe, l|14|l. can now be rewritten as 

R=^((l + 6)y 2 +6(l + 6)xy-l). (34) 

This shows that the boundary R = is given in terms of x and y by (1 + 6)y 2 + 6(1 + 6)xy —1 = and that in region 
I the sign of R must be opposite to the sign of 6 in order to have a positive p. Similarly, in regions II and III, R must 
have the same sign as 6 in order to ensure a positive p. 



10 




FIG. 2: Phase plane of cosmological solutions for uj = and S = 0.1. 



It can be seen that regions II and III are symmetric under a rotation of 7r and a reversal of the direction of the 
trajectories. As region II is exclusively in the semi-circle y ^ all trajectories confined to this region correspond 
to eternally expanding (or expanding and asymptotically static) universes. Similarly, region III is confined to the 
semi-circle y < and so all trajectories confined to this region correspond to eternally contracting (or contracting 
and asymptotically static) universes. Region I, however, spans the y — plane and so can have trajectories which 
correspond to universes with both expanding and contracting phases. In fact, it can be seen from figure □ that, for 
5 = 0.1 all trajectories in region I are initially expanding and eventually contracting. 

It can be seen from figure [5] that in region I the only stable attractors are, at early times, the expanding point 
10 and at late times the contracting point 9. (By 'attractors at early times' we mean the critical points which are 
approached if the trajectories are followed backwards in time). Both of these points correspond to the solution 

a (x rTT^y 

which describes a slow evolution independent of the matter content of the universe. Notably, region I only has stable 
attractor points, at both early and late times, which have been shown to correspond to t — ^constant; region I therefore 
does not have an asymptotic attractor when t — * oo, for the range of 5 being considered. In region II the only stable 
attractors can be seen to be the static point 5 at some early finite time, to, and the expanding matter-driven expansion 
described by point 3 as t — > oo. Conversely, in region III the only stable attractors are the contracting point 4 for 
t — > — oo and the static point 6 for t — > to. 

Figure shows the compactified phase plane for a universe containing pressureless dust and having S = —0.1. 
Figure |3 is split into three separate regions in a similar way to figure [21 with the boundary between the regions again 
corresponding to R = and is given in terms of x and y by (|34|) . Regions II and III again correspond to expanding 
and contacting solutions, respectively. Region I, still has point 10 as the early-time attractor and point 9 as the 
late-time attractor, but now has all trajectories initially contracting and eventually expanding. There are still no 
stable attractors in Region I which correspond to regions where t —* oo. Region II now has point 7 as an early-time 
stable attractor solution and point 1 as a late-time stable attractor solution, corresponding to a — > t * as t — > oo. Point 
3, which was the stable attractor at late times when 6 = 0.1 is now no longer located in Region II and can instead be 
located in region I where it is now a saddle-point in the phase plane. Interestingly, the value of S for which point 3 
ceases to behave as a stable attractor (S = 0) is exactly the same value of S at which the point moves from region II 
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into region I; so as long as point 3 can be located in region I, it is the late-time stable attractor solution and as soon 
as it moves into region I it becomes a saddle-point. At this same value of 5, point 1 ceases to be a saddle-point and 
becomes the late-time stable attractor for region II, so that region II always has a stable late-time attractor where 
t — ► oo. Region III behaves in a similar way to the description given for region II above, under a rotation of 7r and 
with the directions of the trajectories reversed. 

Phase planes diagrams for u> = with values of S other than 0.1 and —0.1, but still within the range 

V73 + 66w + 9w 2 + 3(1 + lu) V73 + 66w + 9uj 2 - 3(1 + u) 

4(4 + 3w) < K 4(4 + 3w) ' 

look qualitatively similar to those above with some of the attractor properties of the critical points being exchanged 
as they pass each other. Notably, for 6 < — | point 3 returns to region II and once again becomes the stable late-time 
attractor for trajectories in that region. The points that are the stable attractors for any particular value of 5 can be 
read off from the tables in the previous section. 

Universes filled with perfect fluid black-body radiation also retain qualitatively similar phase-plane diagrams to the 
ones above; with the notable difference that the point 3 is always located in region II and is always the late-time 
stable attractor of that region. This can be seen directly from the Ricci scalar for the solution Ijl5|l which is given by 

35(1 + 6) 
R - t* 

and can be seen to have the same sign as 5, for S > — 1, and so is always found in region II. 

For a spatially-flat, expanding FRW universe containing black-body radiation we therefore have that l|15|) is the 
generic attractor as t — > co. Similarly, for a spatially-flat, expanding, matter-dominated FRW universe (|15f) is the 
attractor solution as t — > oo; except when — j < S < 0, in which case it is point 1 (aa(5). 

If we require a stable period of matter domination, during which a(t) ~ ti, we therefore have the theoretical 
constraint 5 > (or 6 < — 4). Such a period is necessary for structure to form through gravitational collapse in the 
post-recombination era of the universe's expansion. 

The effect of a non-zero curvature, 0, on the cosmological dynamics is similar to the general relativistic case. 
The role of negative curvature (n = —1) can be deduced by noting that its effect is similar to that of a fluid with 
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to = —1/3. The solution (|15|l is unstable to any perturbation away from flatness and will diverge away from k — as 
t — > oo. This is usually referred to as the 'flatness problem' and can be seen to exist in this theory from the analysis 
of Carloni et. al. [lflj . 



B. Physical consequences 



The modified cosmological dynamics discussed in the last section lead to different predictions for the outcomes 
of physical processes, such as primordial nucleosynthesis and CMB formation, compared to the standard general- 
relativistic model. The relevant modifications to these physical processes, and the bounds that they can impose upon 
the theory, will be discussed in this section. We will use the solutions l|15fl as they have been shown to be the generic 
attractors as t — > oo (except for the case — i < 6 < when u> = 0, which has been excluded as physically unrealistic 
on the grounds of structure formation). 



1. Primordial nucleosynthesis 

We find that the temperature-time adiabat during radiation domination for the solution (|15J) is given by the exact 
relation 



where, as usual (with units K = c = 1 = &b), 



t 2(1+d) = ^ (35) 
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where g is the total number of relativistic spin states at temperature T. The constant A can be determined from the 
generalised Friedmann equation (|13f) and is dependent on the present day value of the Ricci scalar, through equation 
(|10fl . (This dependence is analogous to the dependence of scalar-tensor theories on the evolution of the non-minimally 
coupled scalar, as may be expected from the relationship between these theories H^). As a first approximation we 
assume the universe to have been matter dominated throughout its later history; this allows us to write 

45(l-2<5)(l-25-5<5 2 )\ /2(l + 5) x 2> 



A [ 32(1 - S)g^G ) {-^r) (36) 

where Hq is the value of Hubble's constant today and we have used the solution l|15fl to model the evolution of a(t). 
Adding a recent period of accelerated expansion will refine the constant A, but in the interests of brevity we exclude 
this from the current analysis. 

As usual, the weak-interaction time is given by 

1 

twk 7^- 

The freeze-out temperature, Tj, for neutron-proton kinetic equilibrium is then defined by 

t(T f )=t wk (T f ), 

hence the freeze-out temperature in this theory, with 5 =/= 0, is related to that in the general-relativistic case with 
5 = 0,Tf R ,bj 

„ „ 3(1 + 5) 

Tf = C(Tf R )T3+m (37) 



5(1 + 5) 



where 

r _( (1 ~ <*) \ ( 45 \ ^57T / 3 ff \ (3 + 55 

^(l-2<5)(l-2(5-5<5 2 )y \Z2giPG) \2{l + 5)) ' " 

The neutron-proton ratio, n/p, is now determined at temperature T when the equilibrium holds by 

n ( Am 



- = cxp 
P V T 



13 



where Am is the neutron-proton mass difference. Hence the neutron-proton ratio at freeze-out in the R 1+s early 
universe is given by 

n I Am \ 

exp 



P ' \ c{Tf R y-^ 

where 

25 



3 + 5(5' 

The frozen-out n/p ratio in the R 1+s theory is given by a power of its value in the general relativistic case, 
(?L«V7,by ' 



n In 
p \p 



C(Tf R ) 



GR 



It is seen that when C(T^ R ) £ > 1 (6 < 0) there is a smaller frozen-out neutron-proton ratio that in the general- 
relativistic case and consequently a lower final helium-4 abundance than in the standard general-relativistic early 
universe containing the same number of relativistic spin states. This happens because the freeze-out temperature is 
lower than in general relativity. The neutrons remain in equilibrium to a lower temperature and their slightly higher 
mass shifts the number balance more towards the protons the longer they are in equilibrium. Note that a reduction 
in the helium-4 abundance compared to the standard model of general relativity is both astr ophy sically interesting 
and difficult to achieve (all other variants like extra particl e sp eci es pll 1221 ] , anisotropies [2^, l24l l25l l2q | , magnetic 
fields [13,|2^|, gravitational waves or varying G jlSt l30t l3l| . lead to an increase in the expansion rate and 

in the final helium-4 abundance). Conversely, when C(T^ R ) £ < 1 (5 > 0) freeze-out occurs at a higher temperature 
than in general relativity and a higher final helium-4 abundance fraction results. The final helium-4 mass fraction Y 
is well approximated by 

(1 + n/p) 

It is now possible to constrain the value of S using observational abundances of the light elements. In doing this we 
will use the results of Carroll and Kaplinghat [33] who consider nucleosynthesis with a Hubble constant parametrised 
by 

H ^=(lMevT H - 



Our theory can be cast into this form by substituting 



2 



and 



(1+6) 1 2 

Hi = 1 I J - A-W+sj(lMeV)T^J, 



2 

so, taking g = 43/8, G = 6.72 x 10" 45 AfeV"" 2 and H = 1.51 x 10~ 39 MeV [H, this can bc rewritten as 

(1 + 6)/ 7.96 x 1(T 43 (1 - 5) /2.23 x 1(T 39 \ ^ Tr 

^ = ^( (1-2^(1-2^) ) ( (1 + 5) ) MeV. 

Carroll and Kaplinghat use the observational abundances inferred by Olive et. al. [34[ 

0.228 Y P < 0.248 

2 < 10 5 x— sC 5 
H 

l^HPx— ^3 
H 
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to impose the constraint 

\MeV J 

where H c = 2.6 ± 0.9 x IQ-^MeV at T c = 0.2MeV for 0.5 7710 < 50, or H c = 2.0 ± 0.3 x 10~ 23 for 1 sC -q w ^ 10 
and rji Q is 10 10 times the baryon to photon ratio. 

These results can now be used to impose upon 8 the constraints 

-0.017 s$ 5 sC 0.0012, 

for 0.5 < 7710 ^ 50, or 

-0.0064 sc 5 ^ 0.0012, (40) 

for 1 sC r/10 s$ 10. 



2. Horizon size at matter-radiation equality 

The horizon size at the epoch of matter-radiation equality is of great observational significance. During radiation 
domination cosmological perturbations on sub-horizon scales are effectively frozen. Once matter domination com- 
mences, however, perturbations on all scales are allowed to grow and structure formation begins. The horizon size at 
matter- radiation equality is therefore frozen into the power spectrum of perturbations and is observable. Calculation 
of the horizon sizes in this theory proceeds in a similar way to that in Brans-Dicke theory [85l | . 

In making an estimate of the horizon sizes in R 1+s theory we will use the generalised Friedmann equation, l|13|l . in 
the form 

H +6H R- 6(1 + *) 3(1 - 25) 1¥ P - (41) 

Again, we assume the form (I15f) to model the evolution of the scale factor during the epoch of matter domination. 
This gives the results 

2(1 + 6) 
( t \ 5 

a(t) = a I — 



Ho 

Pm 

R(t) 



2(1 + 5) 

ZHl (1 - 25)(2 - 3(5 - 8<5 2 ) ag 
16ttG (l-<5)(l+5) 2 ^ 
4(1 + 55 + 45 2 ) 
3^2 



during the matter-dominated era. In order to simplify matters we assume the above solutions to hold exactly from 
the time of matter-radiation equality up until the present (neglecting the small residual radiation effects and the very 
late time acceleration). Substituting them into H41[) along with p eq — 2p meq at equality we can then solve for H eq to 
first order in 5 to find 

^hl ~ VVTTI^^ (1 - 2.6865) + 0(6 2 ) (42) 
ao Ho 

where z eq is the redshift at matter radiation equality and H has been treated as an independent parameter. The 
value of 1 + z eq can now be calculated in this theory as 

l + Zeq = ^L. (43) 

PmO 

Taking the present day temperature of the of the microwave background as T = 2.728 ± 0.004K gives 

p r o = 3.37 x 10~ 39 AfeV 4 (44) 
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where three families of light neutrinos have been assumed at a temperature lower than that of the microwave back- 
ground by a factor (4/11)3. Using the same values for G and Hq as above we than find from the above expression 
for p m that 

^-^xl^ C-^P-y ^. (45, 

Substituting (|44|l . I|45|) and (|43J) into (|42J) then gives the expression for the horizon size at equality, to first order in S, 

as 

^tlhl- 155(1 _ 195) +0(6 2 ). (46) 
ao Ho 

This expression shows that the horizon size at matter-radiation equality will be shifted by ~ 1% for a value of 
5 ~ 0.0005. This shift in horizon size should be observable in a shift of the peak of the power spectrum of perturbations, 
compared to its position in the standard general relativistic cosmology. Microwave background observations, therefore, 
allow a potentially tight bound to be derived on the value of <5. This effect is analogous to the shift of power-spectrum 
peaks in Brans-Dicke theory (see e.g. j3^, 

A full analysis of the spectrum of perturbations in this theory requires a knowledge of the evolution of linearised 
perturbations as well as a marginalization over other parameters which can mimic this effect (e.g. baryon density). 
Such a study is beyond the scope of the present work. 

IV. STATIC AND SPHERICALLY-SYMMETRIC SOLUTIONS 

In order to test the R n gravity theory in the weak-field limit by means of the standard solar-system tests of general 
relativity we need to find the analogue of the Schwarzschild metric in this generalised theory and use it to describe 
the gravitational field of the Sun. In the absence of any matter the field equations @ can be written as 

/ R cd R c R d \ ( 1 (1 + 26) \ 

Rab = 6 [ ~R~ S > fl2 ' J \ 9ac9bd + 2 (i-S) 9ab9cd ) ' ^ 47 ' 

We find that an exact static spherically symmetric solution of these field equations is given in Schwarzschild coordinates 
by the line-element 

dr 2 

ds 2 = -A(r)dt 2 + — — + r 2 (d9 2 + sin 2 0d(j> 2 ) (48) 



B(r) 



where 



. , , 25 {1Jr2S) C 

Mr) = r a-*) + - 1T ^ V 



r a 3 *) 
\2 



B(r) = 



(i-sy i , c 



{1-25 + A5 2 )(l - 26(1 + S))\ r {1 -^ s2) 



and C — constant. This solution is conformally related to the Q = limit of the one found by Chan, Home and 
Mann for a static spherically-symmetric space-time containing a scalar-field in a Liouville potential [38|]. It reduces 
to Schwarzschild in the limit of general relativity: 5 — > 0. 

In order to evaluate whether or not this solution is physically relevant we will proceed as follows. A dynamical 
systems approach will be used to establish the asymptotic attractor solutions of the field equations 14711 . The field 
equations will then be perturbed around these asymptotic attractor solutions and solved to first order in the perturbed 
quantities. This linearised solution will then be treated as the physically relevant static and spherically-symmetric 
weak-field limit of the field equations l|47l) and compared with the exact solution l|48|l . 

A. Dynamical system 

The dynamical systems approach has already been applied to a situation of this kind by Mignemi and Wiltshire 
|39|. We present a brief summary of the relevant points of their work in the above notation; for a comprehensive 
analysis the reader is referred to their paper. 
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Taking the value of sign(ii) from (|48|l as sign(— 6(1 + S)/(l —26(1 + 6))) and making a suitable choice of flo allows 
the scalar-field potential 10} to be written as 



35 2 / WirG(l-8) ,\ 

V ® = ~ 8^(1-26(1 + 6)) ^ (V —6* J • (49) 

In four dimensions Migncmi and Wiltshire's choice of line element corresponds to 

ds 2 = e 2U ^ (-dt 2 + f 4 (£K 2 ) + f 2 (£)(dd 2 + sin 2 0d(j> 2 ) (50) 

which, after some manipulation, gives the field equations as 

r „_ 2c\(l - S) 2 + 65 2 v ' 2 - 2A5 2 jC' - 2(1 -26- 86 2 )(' 2 2( 



1-25 + A5 2 
) 2 + 36 2 t 1 ' 2 

Z5 2 (l - 25 + A5 2 ) ' 35 2 



„ _ [1-25- 85 2 )(c\(l - 5) 2 + 35 2 j 2 - 125 2 jC -{1 - 25 - 86 2 )(> 2 ) (1-26(1 + 6)) 2( - 



and 



where 



^ = ~ 3^T^tS) (c?(1 ~ 6)2 + ^ 125Vc ' ~ (1 ~ 25 ~ 862)c ' 2 + (1 ~ 25 + A52)eK) (53) 



C(0 = ^(0+logr(0 

= - (1 ~ 2 3g + g) W ) + 2 logf(0 - + constant. 

Primes denote differentiation with respect to £ and ci is a constant of integration. 
Defining the variables X, Y and Z by 

X = (' Y = r) Z = e c 

the field equations l|51|l and [1 ■"> 2 fl can then be written as the following set of first-order autonomous differential equations 

2c\(l-5) 2 + %5 2 Y 2 -2A5 2 XY - 2(1 - 25 - 85 2 )X 2 2 

X 1 - 26 + 46* Z (M) 

, (l - 26 - 8£ 2 )(c 2 (l - £) 2 + 35 2 Y 2 - 125 2 XY -(1 - 25 - 85 2 )X 2 ) (1 - 26(1 +6)) 2 

35 2 (l - 25 + 4<5 2 ) + 35 2 { ' 

Z' = XZ. (56) 

(The reader should note the different definition of Z here to that of Mignemi and Wiltshire) . As identified by Mignemi 
and Wiltshire, the only critical points at finite values of X, Y and Z are in the plane Z = along the curve defined 
by 

c\(l - 5f + Z5 2 Y 2 - 125 2 XY -(1 - 25 - 85 2 )X 2 = 0. 

These curves are shown as bold lines in figure 0] together with some sample trajectories from equations l|53|land lfS3|l . 
From the definition above we see that the condition Z = is equivalent to fe u = 0. Whilst we do not consider 
trajectories confined to this plane to be physically relevant we do consider the plot to be instructive as it gives a 
picture of the behaviour of trajectories close to this surface and displays the attractive or repulsive behaviour of the 
critical points, which can be the end points for trajectories which could be considered as physically meaningful. The 
dotted line in figure 01 corresponds to the line Y — 2X and separates two different types of critical points. The critical 
points with Y > 2X can be seen to be repulsive to the trajectories in the Z = plane and correspond to the limit 
£ — > — oo. Conversely, the points with Y < 2X are attractive and correspond to the limit £ — > oo. As Z = fe u , all 
critical points of this type in the Z = plane correspond either to naked singularities, f — > 0, or regular horizons, 
f — >constant. 

The two bold lines in figure 0] are the points at which the surface defined by 

c\(l - 5) 2 + Z5 2 Y 2 - 125 2 XY -(1 - 25 - 85 2 )X 2 + (1 - 25 + A5 2 )Z 2 = 
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FIG. 4: The Z — plane of the phase space defined by X ,Y and Z for 5 — 0.1 and ci = 0.5. The bold lines show the critical 
points in this plane and the diagonal lines show the unphysical trajectories confined to this plane. The dotted line is Y — IX 
and separates the critical points where £ — » oo from the points where £ — > — oo 



crosses the Z = plane. This surface splits the phase space into three separate regions between which trajectories 
cannot move. These regions are labelled /, II and III in figure 0] It can be seen from (|53l) that trajectories are 
confined to either regions I or II, for the potential defined by H49fl . If we had chosen the opposite value of sign(iZ) 
in then trajectories would be confined to region III. We will show, however, that region III does not contain 
solutions with asymptotic regions in which f — > oo and so is of limited interest for our purposes. 

In order to find the remaining critical points it is necessary to analyse the sphere at infinity. This can be done by 
making the transformation 

X = p sin 9 cos Y = p sin 9 sin Z = p cos 9 
and taking the limit p — > oo. The set of equations f5K|l and f56")l then give 

f - 24 ^ (1 CQS 2 ^ + 4(5) ( 6 ' 52 cos # 3 " 3S ( 2 ~ 9S ) + C 1 - 5 ( 2 + US )) cos 2d ) 

- (3 - 35(4 - 5(15 - 225 - 32<5 2 )) + (5 - (5(20 - <5(3 + M5 + 32<5 2 ))) cos 29) sin </> 

- 2(18(5 2 (1 - (5(2 + 75)) cos 3cf> - (1 - (5(4 + (5(9 - 265 + 32(5 2 ))) sin 30) sin 2 0) 

and 



deb 



-(6(5 2 (1 - (5(2 + 41(5) - 5(1 - <5(2 + 58)) cos 20)cosec6» sin 



dr 24<5 2 (1 - 2(5 + 4<5 2 ) 

+ 2((1 - (5(4 + (5(9 - 26(5 + 32(5 2 ))) cos 30 + 18<5 2 (1 - (5(2 + 75)) sin 30) sin 9 

- 2 cos 0(4(1 - 25(2 - (5(3 - 28- 4(5 2 )))cosec<9 - (7 - 5(28 + 5(15 - 2(5(43 + 1285)))) sin0)) 

where dr = pd£. These equations can be used to plot the positions of critical points and trajectories on the sphere at 
infinity. The result of this is shown in figure Once again, these trajectories do not correspond to physical solutions 
in the phase space but are illustrative of trajectories at large distances and help to show the attractive or repulsive 
nature of the critical points. The surface at infinity has eight critical points, labelled A-H in figure |SJ Points A and 
B are the end-points of the trajectory that goes through the origin in figure 0] and are located at 

= and 0i (21 = cos^ 1 [ — ^ | (+7r) 

2 A ) V VI - 45 - 12<5 2 + 325 3 + 1005 4 / 
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or, in terms of the original functions in the metric l|5U|) . 

3i 2 3S 2 

f (£ - £1) ci-m-b*') and e u -> - 6) (1_2,5 " 8 ' 52) , 

where ^ is a constant of integration. The points A and B therefore both correspond to £ — » £1 and hence to f — > 0. 

Points C, D, E and _F are the four end points of the two curves in figure 0] and therefore correspond to £ — ► 00 or 
—00 and r — > or constant. 

The remaining points, G and -ff, are located at 



7T 



, 1 _! / 1-25+10S 2 

1.W = 4I+*) and fl=-cos f - § - - - ^ 2 



(1-^+4 ^) /(l — 25 — 2<5 2 ) 1 , { i-28+^)u / (1 — 25 — 25 2 ) 1 , % 

f—a=sr- _ ± ^L___^__ md .-^^ L_^_ , (57) 

where £ 2 is an integration constant, the positive branch corresponds to point H and the negative branch to point 
G. These points are, therefore, the asymptotic limit of the exact solution (|48|) and correspond to £ — > £2 and hence 
f — > 00. 

Whilst it may initially appear that trajectories are repelled from the point H , this is only the case in terms of the 
coordinate £. In terms of the more physically relevant quantity f, the point H is an attractor. This can be seen from 
the first equation in (|57l) . Taking the positive branch here it can be seen that f increases as £ decreases. So, in terms 
of f the points G and H are both attractors, as f — > 00. 

We can now see that in region I all trajectories appear to start at critical points corresponding to either f — ► 
or constant and end at point G where f — > 00. Region 7T appears to share the same features as region J with all 
trajectories starting at either f — > or constant and ending at H where f 00. Region /// has no critical points 
corresponding to f — > 00 and so all trajectories both begin and end on points corresponding to f — > or constant. 

Therefore the only solutions with an asymptotic region in which f — > 00 exist in regions I and II where the potential 
can be described by equation l|49l) . Furthermore, all trajectories in regions I and II appear to be attracted to the 
solution 

ds> = -fAtf + fl-M + 4^(l-2i-2J') dpa + + s . n2 

(1 -<5) 4 
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which is the asymptotic behaviour of the solution found by Chan, Home and Mann [33. We therefore conclude that 
all solutions with an asymptotic region in which f — > 00 are attracted towards the solution (|58|) as f — * 00. 

Rescaling the metric back to the original conformal frame we therefore conclude that the generic asymptotic attractor 
solution to the field equations, 14711 . is 

ds * = + (1-2* + 4^(1-2^-2^ + + s . n2 

(l -sy 

as r — > 00. ft reduces to Minkowski in the 5 — > limit of general relativity. 



B. Linearised solution 

We now proceed to find the general solution, to first order in perturbations, around the background described by 
(|59ll . Writing the perturbed line-element as 

ds 2 = _ r 2S^ {l + v{r))dt 2 + (l- 25 + 45 2 )(l- 25 - 25 2 ) + w{r))dr2 + ^2 + gin 2 (6Q) 

(1 - sy 

and making no assumptions about the order of R the field equations (|47|l become, up to first order in V and W, 

5(l + 25)(l + 25 2 ) (1 + 25 2 ) V 5(1 - 28) W 

(l-6) 2 r 2 + (1-6) T~ 2(1-5) ~ + ~ 

_ 5(1 + 25) R /2 5(1 + 25) R" 35 R' 6(1 + 25) (2 + 5) R' 5(1 + 26) R' 

2 R 2 2(1-5) R 4(1 - 5) R 2(1 - 5) 2 r R + 4(1 - 5) R ' [ ' 

5(l + 28)(l- 25 - 25 2 ) 5(1 + 25) V (2 — 5 + 25 2 ) W V" 

(1 - 5) 2 r 2 (1-5) V + 2(1-6) r 2~ 

36 R' 2 36 R" 5(1 + 25)(2 - 5 + 25 2 ) R' 5(1 + 25) i?' 35 R' 

= 1 1 — — 1 — - — V W (62) 

2 R? 2(1 -5) R 2(1 - 5) 2 r R 4(1-5) R 4(1 - 6) R V ' 

and 

25(3 - 15 + 25 2 + 85 3 ) 2(1 - 26 + 45 2 )(1 - 26 - 26 2 ) TTr V W 

H — — -W 1 

(1 — 5) 2 r 2 (1 — 5) 2 r 2 r r 

m 1 ^ R ' 2 1 5(1 + 25) R" 5(A- 5 + 25 2 + ^) R' 5(1 + 25) R\ r , 6(1 + 25) R' 
Expanding R to first order in V and W gives 



where 

2(1 + 5 + 5 2 ) W 2(1 - 5)(1 + 25 2 ) V 



Itl (l- 25 - 25 2 )r 2 (l- 25 - 25 2 )(l-25 + A5 2 ) r 

(l- 5)(2 - 5 + 25 2 ) W (1-5) 2 yll m) 



(1 - 25 - 25 2 )(1 - 25 + A5 2 ) r (1 - 25 - 25 2 )(1 - 25 + 45 2 ) 
Substituting 1|64|) into the field equations l|61|) . 1|62|) and <|63[l and eliminating R\ using i|65[l leaves 

(l + 5 + 5 2 )(5 - 125+ 125 2 + 45 3 ) W , (16 - 475 + 765 2 - 345 3 - 165 4 + 325 5 ) 



3(l-5) 2 (l + 5) r 6(l-5 2 )(l- 25 + 45 2 ) 

(1 + 5 + 75 2 - 195 3 + 445 4 + 205 5 ) Y (8 - 155 + 185 2 + 165 3 ) , 
3(l-5 2 )(l- 25 + 45 2 ) 7 ~ 6(1 + 5)(1 - 25 + 45 2 ) 



(l- 25 - 25 2 )(5 - 125+ 125 2 + 45 3 ) V (1 - 25 - 25 2 ) 

12(l-5 2 )(l + 5) 7 4(1 -5 2 ) ^ ' 
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(1 + 2S)(1 + 5 + 5 2 )(3 - 45 + 45 2 ) W _ (1 + 2(5) (2 - S + 25 2 )(3 - 45 + 45 2 ) 
3(1 -(S) 2 (l + <5) 7 6(1 -5 2 )(1 — 2c5 + 4<5 2 ) 

(3 - 25+ 175 2 - 45 3 +405 2 ) y (6 - 5 + 25 2 + 205 3 ) 
+ 3(l-5 2 )(l- 25 + 45 2 ) 7 + 6(l + 5)(l- 25 + 45 2 ) 

_ (1 + 2<5)(1 - 25 - 25 2 )(3 - 45 + 45 2 ) ^ (1 + 25)(1 - 25 - 25 2 ) 



12(l-5) 2 (l + 5) r 12(1 -(5 2 ) 

and 

2(8 - 8(5 + 3(5 2 + 10(5 3 - 28<5 4 - 12(5 5 ) W (13 - 225 + 12<S 2 + 26<5 3 - 56<5 4 ) 



3(l-5 2 )(l + 5)(l- 25 - 25 2 ) r 3(1 - <5 2 )(1 - 25 - 25 2 )(1 - 25 + 45 2 ) 



W 



2(4 + 95 2 + 85 3 - 12(5 4 ) Y , (5 - 45 - 45 2 + 125 3 ) , 



3(1 + (5)(1 - 25 - 25 2 )(1 - 25 + 45 2 ) r 3(1 + <5)(1 - 25 - 2<5 2 )(1 - 25 + 45 2 ) 

_ (5 - 45 - 45 2 + 125 3 ) V> (1 + 25) , 
~ 6(l-5) 2 (l + 5) 7 + 6(1 -5 2 y 

where Y = rV and ip = r 3 R[, subject to the constraint 1)650. 

For — < ' 7+ 2o / ^"' > < 5 < — ^ ^he S enera l solution to this first order set of coupled equations is given, in terms 

of \ mid \\ . by 

V(r) = aVi{r) + c 2 V 2 {r) + c 3 V 3 (r) + constant (66) 
W(r) = -aVi(r) + c 2 W 2 (r) + c 3 W 3 {r) (67) 

where 

(1-25 + 45 2 ) 



Vi = 


—r 


(1-5) 




v 2 = 


(1 


+ 25)r 


(1-25 + 45 2 ) 
2(1-5) 


2(2- 


-35 + 


125 2 + 165 3 ) 


w 2 = 


(i 


-25+45 2 ) 


sin(Alogr) 


r 


2(1-5) 


v 3 = 


(1 


+ 25)r 


(1-25 + 45 2 ) 
2(1-5) 


2(2- 


-35 + 


125 2 + 165 3 ) 


w 3 = 


(i 


-25+45 2 ) 


cos(Alogr) 


r 


2(1-5) 



((1 + 25) 2 sin(Alogr) + 2A(1 - 6) cos(Alogr)) 



((1 + 25) 2 cos(Alogr) - 2A(1 - 5) sin(Alogr)) 



and 

V7 - 285 + 365 2 - 165 3 -"805 4 " 

- 2(T~5) • 

The extra constant in (|66|l is from the integration of Y and can be trivially absorbed into the definition of the time 
coordinate. The above solution satisfies the constraint (|65H without imposing any conditions upon the arbitrary 
constants c±, c 2 and c 3 . 

It can be seen by direct comparison that the constant C\ is linearly related to the constant C in l)48|) by a factor 
that is a function of 5 only. The constants c 2 and c 3 correspond to two new oscillating modes. 



C. Physical consequences 



In order to calculate the classical tests of metric theories of gravity (i.e. bending and time-delay of light rays and 
the perihelion precession of Mercury) we require the static and spherically symmetric solution to the field equations 
(J2J. Due to the complicated form of these equations we are unable to find the general solution; instead we propose to 
use the first-order solution around the generic attractor as r — > oo. This method should be applicable to gravitational 
experiments performed in the solar system as the gravitational field in this region can be considered weak and we will 
be considering experiments performed at large r (in terms of the Schwarzschild radius of the massive objects in the 
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system). To this end we will use the solution found at the end of the previous subsection. We choose to arbitrarily 
set the constants C2 and C3 to zero - this removes the oscillatory parts of the solution, and hence ensures that the 
gravitational force is always attractive. This considerable simplification of the solution also allows a straightforward 
calculation of both null and timelike geodesies which can be used to compute the outcomes of the classical tests in 
this theory. 



1. Solution in isotropic coordinates 



Having removed the oscillatory parts of the solution we are left with the part corresponding to the exact solution 
(|48p. Making the coordinate transformation 



C 



(1-25 + 45 2 ) 



r (l-25+4^)/(l-d) = 1 _ ^ f\l (1-25-252) 

I / (1-25 + 452) I 

\ 4^V (1-25-2^) / 

the solution Ij48(l can be transformed into the isotropic coordinate system 

ds 2 = -A{f)dt 2 + B(f){df 2 + f 2 {dd 2 + sin 2 9d<p 2 )) (68) 

where 



2 



23(1 + 2,5) / C \ I C 

A(f) = f \A1-2<S-2S 2 )(1_2S+4,S2) ^ _| - 1 



2(1+45) 
' (1-25+432) 



(1-25+45 2 ) I \ / (i_25 + 45 2 ) 

^-y (1-25-252) J y (i_ 2 5-252) 



and 



B(f) = f V(l-25-25 2 )(l-25+452) M . 



4(1-5) 
(1-25 + 45 2 ) 



Af 

which is. to linear order in C, 



(1-25 + 452) 
(1-25-25 2 ) 



2S(1 + 2S) 2 I (i-<j)(i-2<y) c 

A(f) — f V / (l- 2 «-23 ;i )(l-25 + 452 ) ^ I V M £ - 

1 (1-25 + 4S 2 ) I (1-25 + 452) 



(1-25-252) 



and 



_B(f ) = r V(l-25-252)(l-25 + 452) M 



(1-25 + 4J 2 ) / (1-25 + 452) 
' (1-25-252 ) 



2. Newtonian limit 



We first investigate the Newtonian limit of the geodesic equation in order to set the constant C in the solution l|68|) 
above. As usual, we have 

<j> — 

— 1 oo 

where $ is the Newtonian gravitational potential. Substituting in the isotropic metric 16811 this gives 



1-2 



1-26-25 2 , 3(1-25) 



S(l + 26)f "V (l-<5)(l-8<5 + 4<5 2 )C*^ ^(1-25-252,(1-25+452) ^ 



,/(!- 2<5 -2<5 2 )(1 -2J + 45 2 ) 2^/(1 - 2<S - 2<5 2 )(1 - 26 + 4<5 2 ) s 
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The second term in the expression goes as ~ r 2 +°( s ) and so corresponds to the Newtonian part of the gravitational 
force. The first term, however, goes as ~ f-i+o(<5-) anc j nag n0 Newtonian counterpart. In order for the Newtonian 
part to dominate over the non-Newtonian part we must impose upon 8 the requirement that it is at most 



8-0 



If 8 were larger than this then the non-Newtonian part of the potential would dominate over the Newtonian part, 
which is clearly unacceptable at scales over which the Newtonian potential has been measured and shown to be 
accurate. 

This requirement upon the order of magnitude of 8 allows (|69f) to be written 

where expansions in C have been carried out separately in the coefficients and the powers of f of the two terms. 
Comparison of (|71|l with the Newtonian force law 

_,. Gm 
= — 

allows the value of C to be read off as 

C = -2Gm + 0(8). 



3. Post- Newtonian limit 

We now wish to calculate, to post-Newtonian order, the equations of motion for test particles in the metric ||BSJ|. 
The geodesic equation can be written in its usual form 

dX 2 13 dX dX ~ ' 

where A can be taken as proper time for a timelike geodesic or as an affine parameter for a null geodesic. In terms of 
coordinate time this can be written 



~dt 2 ' 



^ dx 1 * \ dx l dx j 



\\\ also have the integral 



dx 1 dx J 



(72) 



(73) 



where S = — 1 for particles and for photons. 

Substituting H68|) into (|72|l and \T6\ gives, to the relevant order, the equations of motion 



fx 
~dt 2 ~ 



Gm 



dx 



u m 
e r + 4 — e r 



Gm dx dx 
~ er '~dt~dt 

8 



dx 
~dl 



A 52 

4 — e r 

r 



.Gm 



e r + 0(G d m J ) (74) 



and 



dx 
~dt 



= 1-4- 



Gm 



— _ 23 Gm 

j,2S j,l-\-25 



Q(G 2 m 2 ). 



(75) 



(In the interests of concision we have excluded the 0{8 2 ) terms from the powers of r, the reader should regard them as 
being there implicitly). The first three terms in equation (|74ll are identical to their general-relativistic counterparts. 
The next two terms are completely new and have no counterparts in general relativity. The last term in equation 
(I74|l can be removed by rescaling the mass term by m — > m{\ + 8); this has no effect on the Newtonian limit of the 
geodesic equation as any term Gm8 is of post-Newtonian order. 
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4- The bending of light and time delay of radio signals 

From equation (|75|l it can be seen that the solution for null geodesies, to zeroth order, is a straight line that can be 
parametrised by 

x = n(t - t ) 

where n • n = 1. Considering a small departure from the zeroth order solution we can write 

x = n(f - t ) + xi 

where xi is small. To first order, the equations of motion 1)74(1 and ((75(1 then become 

d 2 K „Gm .Gm , , 

— = -2— e r + 4— (n-e r )n (76) 
at r z r z 

and 

dx n Gm 

n - — = -2 . (77) 

dt r K ' 

Equations l(7ti|) and (|77|l can be seen to be identical to the first-order equations of motions for photons in general 
relativity. We therefore conclude that any observations involving the motion of photons in a stationary and spherically 
symmetric weak field situation cannot tell any difference between general relativity and this R 1+s theory, to first 
post-Newtonian order. This includes the classical light bending and time delay tests which should measure the 
post-Newtonian parameter 7 to be one in this theory, as in general relativity. 



5. Perihelion precession 

In calculating the perihelion precession of a test particle in the geometry (|68H it is convenient to use the standard 
procedures for computing the perturbations of orbital elements (see 01 and In the notation of Robertson and 

Noonan |4l| the measured rate of change of the perihelion in geocentric coordinates is given by 

dCo P K J(p + r) . 

— = cos</H -sm<j> (78) 

dt he he 

where p is the semi-latus rectum of the orbit, h is the angular-momentum per unit mass, e is the eccentricity and 1Z 
and J are the components of the acceleration in radial and normal to radial directions in the orbital plane, respectively. 
The radial coordinate, r, is defined by 

(79) 



(1 + ecos^) 

and 4> is the angle measured from the perihelion. We have, as usual, the additional relations 

p = a(l- e 2 ) 

and 

,d(j) 



h = ^Gmp = r 2 -^. (80) 



From 1(74(1 , the components of the acceleration can be read off as 



Gm Gm 2 Gm 2 G 2 m 2 5 6 2 . <5 2 



n = — 5 ^-^ + 4—^ + 4 — = + _^_4_ (81) 



and 



(82) 
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where we now have the radial and normal-to-radial components of the velocity as 

eh . 
v n = — sin 
P 

vj — — (1 + ecos(j>) 
P 

and = v\ + Vj. In writing l|81|) . the last term of l|74|) has been absorbed by a rescaling of m, as mentioned above. 

The expressions l|HT|) and can now be substituted into l(75|) and integrated from cj) = to 2n, using (|79)l and 
(|80ll to write r and dr in terms of </> and d(f>. The perihelion precession per orbit is then given, to post-Newtonian 
accuracy, by the expression 

a(l — e l ) e £ \ Gm ) 

The first term in 1|83|) is clearly the standard general relativistic expression. The second term is new and contributes 
to leading order the term 

2na fl-e 2 " 



Gm V e 



-,2 



8. 



Comparing the prediction (|83[) with observation is a non-trivial matter. The above prediction is the highly idealised 
precession expected for a timelike geodesic in the geometry described by i|68|l . If we assume that the geometry <|68[) is 
a good approximation to the weak field for a static Schwarzschild-like mass then it is not trivial to assume that the 
timclike geodesies used to calculate the rate of perihelion precession ()83|l are the paths that material objects will follow. 
Whilst we are assured from the generalised Bianchi identities 16] of the covariant conservation of energy-momentum, 
T ab b = 0, and hence the geodesic motion of an ideal fluid of pressureless dust, v^u 3 ^ = 0, this does not ensure the 
geodesic motion of extended bodies. This deviation from geodesic motion is known as the Nordvedt effect [42j and, 
whilst being zero for general relativity, is generally non-zero for extended theories of gravity. From the analysis so 
far it is also not clear how orbiting matter and other nearby sources (other than the central mass) will contribute to 
the geometry <|6*51l . 

In order to make a prediction for a physical system such as the solar system, and in the interests of brevity, some 
assumptions must be made. It is firstly assumed that the geometry of space-time in the solar system can be considered, 
to first approximation, as static and spherically symmetric. It is then assumed that this geometry is determined by 
the Sun, which can be treated as a point-like Schwarzschild mass at the origin, and is isolated from the effects of 
matter outside the solar system and from the background cosmology. It is also assumed that the Nordvedt effect is 
negligible and that extended massive bodies, such as planets, follow the same timelike geodesies of the background 
geometry as neutral test particles. 

In comparing with observation it is useful to recast (|83[1 in the form 

. _ 6irGm . 
a(l — e z ) 



where 



A=l+"'<j-fT*. 



This allows for easy comparison with results which have been used to constrain the standard post-Newtonian param- 
eters, for which 

The observational determination of the perihelion precession of Mercury is not clear cut and is subject to a number of 
uncertainties; most notably the quadrupole moment of the Sun (see e.g. 43]). We choose to use the result of Shapiro 
et. al. 

A = 1.003 ± 0.005 (84) 

which for standard values of a, e and m j45j | gives us the constraint 

5 = 2.7 ±4.5 x 10~ 19 . (85) 

In deriving l|84|l the quadrupole moment of the Sun was assumed to correspond to uniform rotation. For more modern 
estimates of the anomalous perihelion advance of Mercury see ^3| . 
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V. CONCLUSIONS 

We have considered here the modification to the gravitational Lagrangian R -> R 1+5 , where 5 is a small rational 
number. By considering the idealised Friedmann-Robertson-Walker cosmology and the static and spherically sym- 
metric weak field situations we have been able to determine suitable solutions to the field equations which we have 
used to make predictions of the consequences of this gravity theory for astrophysical processes. These predictions 
have been compared to observations to derive a number of bounds on the value of 5. 

Firstly, we showed that for a spatially-flat, matter-dominated universe the attractor solution for the scale factor as 
t — > oo is of the form a(t) oc if — j < S < 0. This is unacceptable as sub-horizon scale density perturbations do 
not grow in a universe described by a scale factor of this form. We therefore have the constraints 

6^0 (or 5 < —1/4), (86) 

2(1+5) 

in which case the attractor solution for the scale factor as t — > oo changes to that of the exact solution a(t) oc t 3 

Secondly, we showed that the modified expansion rate during primordial nucleosynthesis alters the predicted abun- 
dances of light elements in the universe. Using the inferred observational abundances of Olive et. al. [34| we were 
able to impose upon 5 the constraints 

-0.017 s$ S sc 0.0012, (87) 

for 0.5 < 7710 ^ 50, or 

-0.0064 sc S ^ 0.0012, (88) 

for 1 7710 10. 

Next, we considered the horizon size at the time of matter-radiation equality. After showing that the horizon size is 
different in this theory to its counterpart in general-relativistic cosmology we discussed the implications for microwave 
background observations. This argument runs in parallel to that of Liddle et. al. for the Brans-Dicke cosmology |3flj . 
The horizon size at matter-radiation equality will be shifted by ~ 1% for a value of 6 ~ 0.0005. 

Finally, we investigated the static and spherically symmetric weak-field situation. We calculated the null and 
timclike geodesies of the space-time to post-Newtonian accuracy. We then showed that null geodesies are, to the 
required accuracy, identical in this theory to those in the Schwarzschild solution of general relativity. The light 
bending and radio time-delay tests should, therefore, yield the same results as in general relativity, to the required 
order. 

Our prediction for the perihelion precession of Mercury gave us our tightest bounds on d. Assuming that Mercury 
follows timelike geodesies of the space-time we used the results of Shapiro et. al. | l ]( to impose upon S the constraint 



6 = 2.7 ±4.5 x 10" 19 . (89) 

This constraint is due to the unusual feature of the static and spherically-symmetric space-time that as r — > 00 it is 
asymptotically attracted to a form that is not Minkowski spacctime, but reduces to Minkowski spacetime as 5 — ► 0. 
Combining the above results we therefore have that S should be constrained to lie within the range 

< 5 < 7.2 x 10" 19 . (90) 



This is a remarkably strong observational constraint upon deviations of this kind from general relativity. 
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the resulting phenomenology of the theory fTBl . Furthermore, the form of the perturbation to general relativity that we 
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